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(X)NVECTIVE AND EQUILIBRIUM RADIATION HEAT-TRANSFER 


PREDICTIONS FOR PROJECT FIRE REENTRY VEHICLE 

By P. Calvin Stairiback 
Langley Research Center 


SUMMARY i . 

C 

Approximate equilibrium flow fields were calculated for the forebody of the 
Project Fire reentry vehicle at zero angle of attack for five points along the 
trajectory. These flow-field calculations were based on wind-tunnel pressure 
and shock-shape data and the assumption of a quadratic variation of pressure and 
velocity across the shock layer. The results of these calculations were ulti- 
mately used to calculate the equilibrium radiation heat -transfer rate to the 
forebody of the vehicle. Convective heating rates were also calculated by using 
the wind-tunnel pressure data. 


INTRODUCTION 


At the present time theoretical estimates for the heat-transfer rates to 
vehicles reentering the atmosphere at very high velocities are subject to rather 
large uncertainties. These uncertainties are particularly true of the component 
of heating due to radiation from the hot gas in the shock layer because of the 
lack of accurate radiation intensity data. One purpose of Project Fire is to 
measure the radiation heating rate (both total and spectral) and total heating 
rates (convective plus radiative minus reradiative) experienced by a body reen- 
tering the atmosphere at a velocity of 37>O00 feet per second. 

Since basic air radiation data can be greatly influenced by foreign parti- 
cle contamination, the vehicle is designed to obtain radiation and total heating 
rates in a "clean" atmosphere. This design is accomplished by utilizing a mul- 
tiple layer construction for the forebody, which will receive the majority of 
the heat load. Three metal (beryllium) calorimeters are used to insure a clean 
atmosphere for a suitable testing time. Two ablating heat shields, which can 
be removed at a selected t im e by an explosive mechanism, are interspaced between 
the calorimeters to protect them until the selected testing time is reached on 
the trajectory. 

With sufficiently extensive and accurate flight data, it should be possible 
to evaluate the various theories for calculating convective heating rates and 
the various radiation intensity data. The purpose of the present report is to 
present approximate equilibrium flow fields and the resultant heat-transfer 
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calculations for the Project Fire reentry vehicle at several points along its 
proposed trajectory. The results of these simplified calculations are compared 
with the results of more exact calculations to determine if the relatively sim- 
ple analysis can provide heating estimates of usable accuracy. 


SYMBOLS 


A 

area 




ai,a 2 ,a 5 

constants 

in 

eq. 

(i) 

bi,b 2 ,b3 

constants 

in 

eq. 

(2) 

c l> c 2 ,c 3 

constants 

in 

eq. 

(2) 


drag coefficient 



h altitude 

h,k,Z coordinates of point on body in X,Y,Z coordinate system 

j specific radiation intensity 

Mi mass flow rate into shock layer 

Mo mass flow rate out of shock layer at 9 

N = n/A 


n 

n 

P 

i 

*t 

q 

ic 


normal distance from body surface 
number of subdivisions in Simpson's rule 
pressure 

stagnation pressure behind normal shock 

heat-transfer rate 
convective heat-transfer rate 


q R 

R 

R c 


radiative heat -transfer rate 
radius 

corner radius 


2 


effective radius 


R eff 

R n nose spherical radius 

r cylindrical radius of body (see fig. 4) 

r mflY maximum cylindrical radius of body 

r distance between point on body and dV 

s surface distance from stagnation point 

T temperature 

t time (t = 0 at 400,000 ft) 

I 

U velocity 

V gas cap volume 

i W vehicle weight 

X,Y,Z body-axis rectangular coordinate system 

i 

x,y,z rectangular coordinate system, origin at point P on body 
i 7 reentry angle 

' A local shock standoff distance from body 

£ * ( P “/ P s)g_o 
£ vorticity 

9 angle measured from center line of body 

9 C angle measured from forebody-corner line of tangency to line passing 

through any point on corner 

l 9p angle measured from center line of body to line passing through 

point P 

0g angle measured from center line of hypothetical sphere 

0ij angle measured from center line of body to forebody-corner line of 

tangency 

p density 

p s i density at sea level 

i 3 



0 




angle measured from body normal to flow-field velocity 
co angle measured from body normal to direction of dV 

Subscripts: 

b at body surface 

s behind shock, or of shock 

t at stagnation point 

oo free stream 


FLOW-FIELD ANALYSIS 


Approximate inviscid flow fields were calculated for the Project Fire 
reentry vehicle at the following trajectory points, which were taken from the 
latest trajectory analysis available at the start of the present investigation: 


Altitude, 

ft 

Velocity, 
ft/ sec 

Time, 

sec 

Configuration 

260,460 

37,058 

15 

1 

218,000 

36,700 

19 

1 

166,000 

34,800 

25 

2 

147,000 

30,100 

27.6 

2 

120,000 

19,500 

32 

3 


The initial conditions for the trajectory from which the points were taken are: 
Uoo = 37,000 ft/sec; h = 400,000 ft; 7 = -15°; and W/C^A = 35 Ib/ft2. The con- 
figuration numbers represent the shapes of the various beryllium calorimeters 
which, with the ablation shields, form the forebody of the vehicle. (See 
fig. 1.) The altitude and velocity at t = 25 seconds are the approximate con 
ditions where the peak total heating rate is expected. A schematic of the heat 
pulse expected for the vehicle and the periods when the beryllium calorimeter 
and radiation sensors are expected to obtain useful data are shown in figure 2. 

The basic assumptions made during the flow-field analysis were: 

(1) Equilibrium flow 

(2) Pressure distribution known (based on wind-tunnel data) 

(3) Shock shape known (based on wind-tunnel data) 
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It can be expected that nonequilibrium effects will be small except possibly at 
the highest altitude noted in the previous table; therefore, equilibrium flow 
was assumed to simplify the analysis. 

Pressure distributions were obtained from low-enthalpy tunnel investiga- 
tions at a Mach number of 8, and in order to apply them to flight conditions, 
it was necessary to assume that the pressure distribution, with respect to the 
stagnation -point value, was invariant with gas conditions. Because of minor 
changes in the final vehicle shape, some adjustment of the pressure distribu- 
tion was required to account for the differences in the models and the final 
vehicle forebody. These adjustments were minor and should have negligible 
effects on the results. The actual pressure distributions used throughout the 
present analysis for the three beryllium calorimeters are presented in figure 3* 

The shocks for configurations 1 and 2, obtained from tunnel investigations 
at a Mach number of 8, were concentric with the forebody to the line of tangency 
between the forebody and the corner. Concentricity of the shock was also found 
to exist for configuration 2 in a high-enthalpy expansion tube. This concen- 
tricity of the shock was assumed for the third calorimeter in lieu of a measured 
shock shape. With this assumption, the shock standoff distance was the only 
unknown required to define the shock envelope over the forebody. The shock 
standoff distance was obtained from the correlation presented in reference 1, 
which can be used to obtain a relationship between and R e ff> and the fol- 
lowing equation for shock standoff distance: A/R e ff «= O.S^p^p^. Downstream 

of the conical surface defined by the forebody-corner tangency line and noted 
by its half -angle 0<p (see fig. k ) , the shock was assumed to be described by a 
second-order curve given as: 


i 

i 




+ Gi? 


This curve was required to pass through the proper point along the line at 


( 1 ) 


0ip 


with the correct slope and to pass through one additional point near the most 
rearward region that could influence the radiation heating up to and including 
the forebody-corner line of tangency. This latter point was obtained from the 
measured shock shape for configuration 2. 

Since the pressure distribution along the body and the shock shape were 
known, the fluid conditions along the body and downstream of the shock could be 
determined. In order to determine the fluid properties in the inviscid shock 
layer, it was assumed that the pressure and velocity distributions across the 
layer along a normal from the surface could be expressed in terms of a second- 
order polynomial of the form: 


*b 



+ b2N + bj 


U 

% 


= c ± N 2 + 


C2N + 


( 2 ) 
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The three unknown coefficients in each of these equations were determined from 
the conditions at the wall, the conditions downstream of the shock, and the 
gradient at the wall. The pressure gradient at the wall, obtained from the 
inviscid momentum equation, is: 



(3) 


The velocity gradient can be obtained from the expression for the vorticitj 
of a fluid behind a shock; this equation for vorticity is given in reference 2 
as: 


(1 - e)2 

Pb T b 

**00 

€ 

( PsTs )e=o 

Rs 2 




The velocity gradient normal to a streamline is: 


dU _ ► U 

dn “ R 


( 5 ] 


Combining equations (4) and ( 5 ) gives the normal velocity gradient along the 
body surface as: 



(1 - e) 2 

’ 0b T b 

Ur 

00 

€ 

( p s T s)q_o 

— -J 

R s 2 


Hb 

*b 


( 6 : 


The values for 




Pb, p s , 1^, and U s 


permit the coefficients 


in equations (2) to be evaluated. With the pressure and velocity known and 
constant total enthalpy assumed throughout the shock layer, the other equilib- 
rium fluid properties can be obtained from suitable gas tables or charts for 
air. (See, for example, ref. 3 .) With the local fluid properties defined, the 
absorption or emission coefficient can be obtained from tabulated data. For 
the present analysis the emission coefficients were obtained from reference 4 

for = 10“^. For — < 10 the emission coefficients were obtained from 

p s l p s l 

reference 5 since values for this low-density range were not included in 
reference 4. 


Since the equations for the pressure and velocity distributions are writte 
in terms of a coordinate system normal to and along the body surface, it is con 
venient to divide the flow-field analysis into two parts, that is, flow over th 
forebody ^9 < 0ip^ and flow around the corner ^9 > 0^. (See fig. 4.) Some 
consideration must be given to the value of R^ since at the line of tangency 


6 


\ 


between the forebody and corner the body radius changes and results in a dis- 
continuous change in the normal pressure and velocity gradients at 9 * 0 T . 
Since the flow is subsonic, it does not appear reasonable from a physical view- 
point to have this discontinuity. Therefore, in order to eliminate this dis- 
continuity in the gradients normal to the body, an effective body radius was 
calculated. This calculation was made by assuming that the effective radius of 
the body is equal to the radius of a sphere which has the same pressure gradi- 
ent along its surface at equal angles with respect to the free-stream velocity. 
A Newtonian pressure distribution was assumed to exist over the hypothetical 
sphere. With these assumptions, the effective body radius becomes: 



2 cos 9c? sin 9c 


where 


( p K) 


a (»At.) 

d(s/R) 


(7) 


is determined from the measured pressure data for the body. 


a(s/R) 

This effective radius was used throughout the analysis for R^ in equations (3) 

and (6) . It should be noted that this definition of R e ff should give the same 
values for R e ff at the stagnation point as given in reference 1. A plot of 
R e ff and R^> for configuration 2 is given in figure 5 to show the difference 
in the two quantities. Figures 6 to 10 present the temperature, density, and 
specific radiation intensity calculated by the previously discussed method for 
both the forebody and corner flow fields. In general, the fluid properties 
vary fairly uniformly over the forebody and around the corner up to a value of 
9 C of about 30° to 40°. For larger values of 9 C the flow-field properties 
appear to be somewhat erratic, and this can probably be attributed to the sim- 
plicity of the present method. 


Except for the gas in the immediate vicinity of the point for which the 
radiation heating rate is being calculated, the irregularity of the flow-field 
fluid properties for large values of 9 C has little influence on the radiation 
heating rate to the forebody due to the distances and view angle involved. The 
limit of the flow field that is viewed by the forebody-corner tangency point is 
noted on the curves in figures 6 to 10 by the dots. 


It should be pointed out that although no iterations were made during this 
analysis, the method employed could be revised to permit iteration on the shock 
location by balancing the mass flow within the shock layer. In order to obtain 
some quantitative evaluation of the self consistency of the calculated flow 
fields, a mass-flow balance was made at 8 * 15° • The results of this balance 
are shown in figure 11 and indicate that based on the criterion of mass flow, 
the flow fields appear to be only fairly accurate. 

The mass-flow balance could be influenced by three things: the calculated 

density and velocity distributions across the shock layer, the shock standoff 
distance, and the flow angularity with respect to the area over which the mass- 
flow balance was calculated. The variation of the flow angle with respect to 
the normal from the body surface was assumed to be linear between the value at 
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the body (0 = 90 °) a- 11 *! "the value calculated behind the shock for an oblique 
shock. The density and velocity distributions across the shock layer depend on 
the value of R e ff • A comparison between the mass-flow balance for R = R e ff 
and R = R n in equations ( 5 ) and ( 6 ) indicates that there is little difference 

between the two methods. The mass -flow balance parameter ^1 - for the 

flow field calculated for R = R e ff was 0.132 and for R = R n was 0.1435. 
Therefore, if the flow-angularity assumption is reasonable, the mass-flow 
deficit must be predominately due to the shock standoff distance which is too 
small. 


RADIATION HEATING CALCULATIONS 


The 
gas that 


quantity of energy received at a point P from a surrounding radiating 
is transparent, non scattering, and in equilibrium can be expressed as 


q.R = 



1 


cos dV 
r 2 


(8) 


Equation ( 8 ) is written for a coordinate system with its origin located at the 
point P receiving the radiant energy. It is convenient to express this equa- 
tion in a body-axis coordinate system; the details of this transformation are 
outlined in the appendix. Applying the resultant equation to a portion of a 
sphere at zero angle of attack gives: 


riy- J Rx - h)cos 9p + (Y - k) sin e^l I dX dY dZ 
^ ' 


[(X - h ) 2 + (Y - k ) 2 + (Z - l ) 2 ] 5 ^ 2 


The limits of integration are as follows: 
for Z, 


= tan 0 ij 

A 




V 

— > tan 9n 






( 9 ) 


( 10 ) 


8 



for Y (upper) , 



for Y (lower), 



(ID 


( 13 ) 


9 



and for X (upper), 



Equation ( 9 ) was integrated by- 
Simpson's male where the number of subdi- 
visions, rT, along the X coordinate axis 
was taken as given in table 1. Along the 
Y- and Z-axes, rT = 100. 

An indication of the convergence of 
the integration to a given value as n 
increases is shown in figure 12 for 
t = 25 seconds; the solid curves represent 
the envelope of the maximum and minimum 
variation of the radiation heating with 
n. If it is assumed that the correct value 
of qp is the average value obtained with n = 100 and n = 102, the value for 

n = 100 is in error only about 1 percent. 

The number of subdivisions, n, along the X-axis was varied from n = 10 
at 8p = 0 to !T = 20 for the largest value of 8p since the absolute value 

of ZsX increased with increasing 8p. 

The radiation heat-transfer distribution and the stagnation-point heating 
rate are presented in figures 13 and l4. All the distributions are similar with 
the radiation heating decreasing to about 0.2 to 0.4 of the stagnation -point 
value at 9 = 8«p . The radiation heat-transfer distribution for the semi- 
infinite slab approximation is also presented in figure 13 • This approximation 
is similar to the slab approximation usually used at the stagnation point except 
that an average value of the specific intensity obtained from the values behind 
the shock and at the wall is used with the local shock standoff distance. The 
results indicate that this approximation gives an excellent indication of the 
distribution when one considers its simplicity. This result is in agreement 
with a more detailed slab approach made in reference 6. 


TABLE 1 


9p, <leg 

n 

0 

10 

5 

12 

10 

14 

15 

16 

20 

18 

>20 

20 
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The mass-flow discrepancy noted previously in the section entitled "Flow- 
Field Analysis" will influence the magnitude of the radiation heating since the 
error is believed to be predominantly due to the shock standoff distance being 
too small by approximately the same percentage as the mass -flow balance. If the 
slab approximation is considered, the absolute magnitude of the heating would 
increase in proportion to the mass-flow error, but the distribution would be 
essentially unchanged. 

The stagnation-point radiative heating rates shown in figure l4 are rela- 
tively low for all points investigated except for t = 25 seconds where the 
heating rate is about 400 Btu/sec-ft^. Also shown in figure l4 is the slab 
approximation for the stagnation-point heating rate. 

It should be noted that the radiative heating rate is dependent on the 
specific-radiation-intensity data used. To date this gas property is not known 
to great accuracy and different values are given by different investigators. 

For the most part the radiation intensity data used for this investigation were 
taken from reference 4 since it was generally accepted that the best data avail- 
able at the start of this investigation are given in this reference. Recently 
it has been suggested that these intensities are low. In order to show the 
influence of the various radiation intensity data, the stagnation-point heating 
rate and the forebody distribution for t — 25 seconds were calculated by 
using the data of references 4, 5> and 7> The results of these calculations 
are presented in figure 15* The stagnation-point heat- transfer rates range 
from 4l0 Btu/sec-ft^, if the data of reference 4 are used, to 946 Btu/sec-ft^, 
for the data from reference 7- There was no change in the radiation heating 
distribution obtained by using the various specific-intensity data. 

The difference between the heating estimates can be partly attributed to 
the larger wavelength range considered for calculating the specific intensities 
in reference 7* The wavelength range for the data of reference 7 was from 0.05 
to 10 microns whereas the range was from 0.l6 to 10 microns for reference 4. 

It is expected that the Project Fire flight data will reduce the uncertainty 
in predicting the radiation heating. It should be noted, however, that most of 
the added energy in the wave length increment from 0.05 to 0.l6 micron con- 
sidered in reference 7 will not be transmitted to the vehicle radiation sensors 
since the quartz windows in the forebody absorb most of the radiation below 
about 0.18 micron. Thus, this increment in radiation heating must be, in some 
way, inferred from convective-heating calculations and measured total heating 
rates. 

In addition to this complication, self absorption of radiant energy within 
the gas cap might become significant, particularly in the small wavelength 
region considered in reference 7* 


CONVECTIVE HEATING 


The convective heat-transfer distribution around the vehicle was calcu- 
lated by using the correlations presented in reference 8. The results of the 
calculations of the heat-transfer distribution are presented in figure 16 . 
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This figure shows that the heating rate is essentially constant over the fore- 
body; the high heating rate sometimes encountered at the outer edge of blunt 
bodies does not occur. The variation at the outer edge of the body is never 
greater than ±10 percent of the stagnation-point value. 

The stagnation -point heat -transfer rate was calculated from the correlated 
results of reference 9 for the high-enthalpy case. The results of these calcu- 
lations are presented in figure 17 • The maximum convective stagnation-point 
heating rate calculated was 7^0 Btu/sec-ft^. This convective heating rate is 
almost twice the radiative heating rate calculated by the data of reference k 
and about 20 percent less than the radiative heating rate calculated by the 
data of reference 7- Thus, it is not clear whether the convective heating rate 
will dominate the peak stagnation -point heating to the vehicle or whether radi- 
ative and convective heating will contribute equally to the total peak heating 
rate. The radiative and total heating rates are also plotted in figure 17- 
The maximum total heating rate to the stagnation point will range from about 
1150 to about I 69 O Btu/sec-ft^ depending upon the radiation data used. 


COMPARISON OF PRESENT RESULTS WITH MORE EXACT ANALYSES 


The present results are compared with the results of more exact calcula- 
tions made by the General Electric Company, Lockheed Missiles and Space 
Company, and Philco Corporation under contract to the NASA Office of Advanced 
Research and Technology. The results of these comparisons are shown in fig- 
ure l 8 . The figure indicates that there is considerable disagreement between 
the four sets of results. This is particularly true for the convective and 
radiative heat-transfer distributions and for the stagnation-point radiative 
heating rate except when the same specific radiation intensity data are used. 

The stagnation-point convective heat-transfer rate of the present report 
compares favorably with the contractor results except for t = 15 seconds where 
the present results are about 50 percent of the highest estimate. The present 
estimates of the stagnation-point radiative heating rate for t ** 25 seconds 
agree well with the other results provided the same specific radiation intensity 
data are used. There is a great deal of difference between the stagnation-point 
radiative heating rate for t = 15 seconds where nonequilibrium effects can 
influence the radiative heating rate. The present results for this case are in 
fair agreement with Philco' s results. 

It is somewhat surprising that the convective and radiation estimates of 
heat-transfer distributions differ by such a large amount. The present results 
fall between the maximum and minimum values predicted by the contractors. The 
present radiation distribution agrees closely with Philco' s results whereas the 
present convective distribution agrees fairly well with General Electric's 
results. From an overall point of view, the simplified analysis of the present 
report gives results which compare favorably with the results obtained from 
more exact calculations. 
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CONCLUSIONS 


Calculations of the approximate equilibrium flow fields for the forebody 
of the Project Fire reentry vehicle at zero angle of attack and the resulting 
convective and equilibrium radiation heat-transfer calculations permit the fol- 
lowing conclusions: 

1. The radiation stagnation-point heating rate is greatly influenced by 
the specific-radiation-intensity data used and ranged from a low of 

4l0 Btu/sec-ft^ to a high of 946 Btu/sec-ft^. 

2. The radiation heat- transfer rate over the forebody decreases fairly 
rapidly from the stagnation point to the outer diameter of the vehicle. The 
heating rate at the forebody- corner tangency point is only about 0.2 to 0.4 of 
the stagnation -point value. 

5. The peak convective stagnation-point heating rate is about 
7^0 Btu/ sec-ft2. This rate combined with the peak radiation heating results 
in a peak total heating rate that ranges from 1150 to 1690 Btu/ sec-ft^ 
depending on the radiation data used. 

4. The convective heat-transfer rate over the forebody is essentially con- 
stant, never varying more than ±10 percent of the stagnation -point value. 

5. For the case considered, the radiative heat-transfer distributions 
obtained by a simple slab approximation are in good agreement with those 
obtained by integrating the specific intensity over the calculated gas cap 
volume . 

6. In general, the simplified analysis of the present report gives results 
which compare favorably with the results obtained from more exact calculations. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., October 21, 1964. 


13 



APPENDIX 


RADIATION HEATING EQUATION 


The heat-transfer rate to a point P from a surrounding radiating gas 
that is transparent, nonscattering, and in equilibrium can be expressed as 


<lR 


-L 


j cos a) dV 
r 2 


(Al) 


This equation is written for a coordinate system with its origin located at the 
point P receiving the radiant energy. It is convenient to express this equa- 
tion in a body-axis coordinate system. First, expressing the equation in terms 
of a rectangular Cartesian coordinate system with the origin at the point P 
gives 


*R 


. rr n{\ * 0 

JJJ / o n 


cbc dy dz 


r2 j. , r 2 ^ .2) 5 / 2 


(A2) 


x. + y + z c 


where : 

(1) x is perpendicular to the body surface at P. 

(2) y and z lie in the tangent plane at P. 

( 5 ) y lies in plane OAFB. (See fig. 19>) 

In terms of a body-axis coordinate system equation (A2) becomes: 




ffl 


j|j(X - h) cos ai + (Y - k)cos P]_ + (Z - l)cos 7-^ | dX dY dZ 


[a(X - h) 2 + B(Y - k) 2 + c(Z - l ) 2 + 2D(X - h)(Y - k) + 2E(X - h)(Z - l) + 2F(Y - k)(Z - l)] 


3/ 2 


(A3) 


where 


pop 

A = cos c O]_ + cos 02 + cos o^ 

B = cos 2 P2 + cos 2 P 2 + cos 2 pj 

C = cos 2 7-^ + cos 2 72 + cos2 ^3 


D 

= COS 

a l 

cos 

Pi 

+ 

cos 

°2 

cos 

P 2 

+ 

cos 

a 5 

cos 

h 

E 

= COS 

a l 

cos 

71 

+ 

cos 

02 

cos 

72 

•f 

cos 

a 3 

cos 

7 3 

F 

= cos 

Pi 

cos 

71 


cos 

P 2 

cos 

72 

•f 

cos 

h 

cos 

7 5 v 




(A4) 


14 



cos ar|_ : cos Pi : cos 7i = 


bf 

M. 


bf 

M. 


b f 
bZ 


/Sf\ 2 + /8f\ 2 + (bf\ 2 

vW/ W/ Vaz/ 




h cos 8i 
cos Op — 

f 


h cos a-\ + l cos 7i > 
cos p 2 = i— i±) 


cos 7 2 = - 


l COS Pi 

"IT” 




(A5) 


(A6) 


cos a? 


cos 


p 5 = 


l (cos 2 Pi + cos 2 7i) + k cos a l cos 7i 

fi 

cos Pi(h cos 7i - l cos ai) 

Ji 






cos 7- 


[h (cos^ai + cos 2 Pi) + i (cos ai cos 7i) 

fi J 


(A7) 


a = (h 2 + l 2 ^cos 2 Pi + cos ai + l cos 7i^‘ 


(A8) 


6 « [l (cos 2 p x + cos 2 7 i) + h^cos ai cos 7 ijj + cos 2 Pi^h cos 7 i - l cos a^ 2 

+ jji^cos^i + cos 2 Pi^ + i cos ai cos 7iJ^ (A9) 

and h, k, and l are the coordinates of the point P in the X,Y,Z coordi- 
nate system. The body surface is given by 


f(X,Y,Z) = 0 


(A10) 
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For the present case, the portion of the body for which the radiation 
heating was calculated was a portion of a sphere and the angle of attack was 
zero. Therefore, h = R^ cos 0p, k * R^ sin 0p, and 1 a 0; and 


COS 

= cos 

0 p 

COS 

ctg = -sin 0p 

COS 

o 

11 

dP 

cos 

= sin 

0 p 

COS 

^2 = cos 0p 

COS 

P3 - 0 

cos 

71-0 


COS 

o 

11 

CV1 

cos 

7 5 = X 


Also, A = B = C = 1 and D = E = F = 0. Thus, the equation for qp becomes: 


j(p,T) |[(X - h)cos 0 P + (Y - k) sin 0 p]|dX dY dZ 
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Predicted total heating rate, Btu/sec-ft* 


Total heating measurement 
Radiative heating measurements 


14 x 10 


t = 25 sec ; radiation 
experiment - no. 2 
data period ends 


t = 24 sec; first 
2 Q- ablative layer 
is j ettisoned 


t = 27 sec; calorimeter 
experiment - no. 2 
data period ends 


g. t = 20 sec; calorim- 
eter experiment - j 
no. 1 data period 
ends / 


t = 31 sec; second 
ablative layer 
is jettisoned 


t = 15 sec; radi-J 
ation experi- J 
ment - no. 1 M 
data period M 
ends \M- 


t = 32 sec ; radiation 
experiment - no. 3 
data period ends 

t = 40 sec; calorimeter 
L experiment - no. 3 
data period ends 


Reentry time, sec 


Figure 2.- Primary data periods showing total and radiative stagnation-point heating 
measurements. (At start of reentry, t = 0: all times are estimated.) 











Temperature, 


106 x 10* 




(a) Forebody temperature. 

Figure 6.- Flow-field fluid properties for configuration 1. t = 15 sec; h = 260,460 ft 

U “ = 37,058 ft/sec. (Dots denote limit of flow field viewed by forebody-corner tan- 
gency point.) 

















, w/cm -steradian 



(f) Comer specific radiation intensity. 
Figure 6.- Concluded. 




(a) Forebody temperature. 


Figure 7 .- Flow-field fluid properties for configuration 1. t = 19 sec; h = 218,000 ft 
Uoo = 36,700 ft/sec. (Dots denote limit of flow field viewed by forebody-corner tan- 
gency point.) 
















(a) Forebody temperature. 


Figure 8.- Flow-field fluid properties for configuration 2. t = 25 sec; h = 
Uco = 5^*800 ft/sec. (Dots denote limit of flow field viewed by forebody • 
gency point.) 


= 166,000 ft: 
‘Corner tan- 









, w/cm -steradian 



(c) Forebody specific radiation intensity. 
Figure 8.- Continued. 
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Temperature, 







(a) Foretiody temperature. 


Figure 9*- Flow-field fluid properties for configuration 2. t = 27.6 sec; h = 147,000 ft 
U ro = 50,100 ft/sec. (Dots denote limit of flow field viewed "by foretody-corner tan- 
gency point. ) 
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Temperature, 


£ 

Eh 



(a) Forebody temperature. 

Figure 10.- Flow-field fluid properties for configuration 3 . t = 32 sec; h = 120,000 ft; 
- 19^500 ft/sec. (Dots denote limit of flow field viewed by forebody-corner tan- 
gency point . ) 
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Figure 11. - Mass-flow balance for shock layer. 



























Figure 14. - Radiation heat transfer to stagnation point. 




Reference 4 


Reference 5 


Reference 7 


(a) Stagnation-point value. 

Figure 15.- Comparison of radiation heating using various values for specific intensity, 
t = 25 sec; h = 166,000 ft; = 54,800 ft/sec. 






(a) t = 15 sec; h = 260,460 ft; u*. = 37,058 ft/sec 
Figure 16.- Convective heat-transfer distribution. 





19 sec; h = 218,000 ft; U*. = 56,700 ft/i 


Figure l6.- Continued. 




25 sec; h = 166,000 ft; U„ = 34,800 ft/ 


Figure 16 .- Continued. 






a) t = 27.6 sec; h = 1^7,000 ft; = 31,100 ft/ 
Figure l6.- Continued. 



32 sec; h = 120,000 ft; U*, = 19,500 ft/ 


Figure 16.- Concluded. 





Convective 

Radiative, reference 4 



Figure 17- - Total heat-transfer rate to stagnation point. 







(a) Stagnation-point convective heating rate. 

Figure l8.- Comparison "between heating rates and heating distributions. 
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(c) Convective heat-transfer distribution. 
Figure 18.- Continued. 







